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Abstract 

. , We establish the global existence and stability of a three-dimensional supersonic conic shock wave for a 

■^-H ' perturbed steady supersonic flow past an infinitely long circular cone with a sharp angle. The flow is described 

(-H , by a 3-D steady potential equation, which is multi-dimensional, quasilinear, and hyperbolic with respect to 

the supersonic direction. Making use of the geometric properties of the pointed shock surface together with 
the Rankine-Hugoniot conditions on the conic shock surface and the boundary condition on the surface of the 
cone, we obtain a global uniform weighted energy estimate for the nonlinear problem by finding an appropriate 
multiplier and establishing a new Hardy-type inequality on the shock surface. Based on this, we prove that a 
multi-dimensional conic shock attached to the vertex of the cone exists globally when the Mach number of the 
^ ' incoming supersonic flow is sufficiently large. Moreover, the asymptotic behavior of the 3-D supersonic conic 

^ : shock solution, that is shown to approach the corresponding background shock solution in the downstream 

domain for the uniform supersonic constant flow past the sharp cone, is also explicitly given. 

I Keywords: Steady potential equation, supersonic flow, multi-dimensional conic shock, global existence, 

■ Hardy-type inequality, tangential vector fields. 
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§1. Introduction 
• 5—1 . 

. In this paper, we are concerned with the multi-dimensional steady and supersonic conic shock wave problem 

^ ' for a perturbed incoming supersonic flow past an infinitely long circular cone. This problem is fundamental in 

gas dynamics, for instance, for the supersonic flight of projectiles and rockets. It is also one of the basic models 
for the discussion of the theory of weak solutions to quasilinear hyperbolic equations and systems in multi- 
dimensions (see [3], [22]-[23], [30], [37]). Under suitable assumptions on the incoming supersonic flow with 
a small spherically symmetric perturbation and a spherically symmetric pointed body or artificial boundary 
conditions on the conic surface, there is an extensive literature studying supersonic flow past a pointed body 
(see [5], [7]-[10], [18], [21], [31], [35]-[36], and the references therein). The first rigorous mathematical analysis 
was given in [7] by Courant and Friedrichs, who proved that, for a uniform supersonic fiow (0, 0, go) with 
constant density po > which approaches from minus infinity, when the fiow hits the sharp circular cone 
\J x\ + x\ = &o2;3, &o > 0, in direction of the xs-axis (see Figure 1 below), then there appears a supersonic 
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conic shock \lx\-\- x\ = SqXz-, sq > ^Oi attached to the tip of the cone provided that 60 is less than some critical 
value 6* > 0, which is determined by the parameters of the incoming flow. When the incoming supersonic flow 
is niulti-diniensionally perturbed, the basic problem of both mathematical and physical relevance that naturally 
arises is whether such a conic shock is globally stable. Or else, do there appear new shocks or other complicated 
singularities in the downstream domain? Here, we will focus on this problem when the Mach number of the 
incoming supersonic flow is appropriately large. It will be shown that a global supersonic conic shock exists 
uniquely in the whole space and that there are no other singularities between the conic shock and the conic 
surface for a multi-dimensionally perturbed supersonic polytropic gas past the sharp cone \/ x\-\- x\ = box^ 
(see Figure 2 below). This result agrees with physical experiment and numerical simulations. 
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Figure 1. A uniform supersonic flow past a sharp circular cone 




Figure 2. A multi-dimensionally perturbed supersonic flow past a sharp circular cone 
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In this paper, we will use the potential equation to describe the motion of a supersonic polytropic gas (this 
model is also favored in [22]-[23], [30]), where polytropic gas means that the pressure P and the density p of 
the gas are related by the equation of state P = Ap"' , with A > being a constant and the adiabatic constant 7 
satisfying 1 < 7 < 3 (for air, 7 ~ 1.4). Let $(a;) be a potential of the velocity u = (mi, ■U2, Ms), i-e., Ui = 
Then it follows from Bernoulli's law that 

i|V$|2 + /i(p) = Co. (1.1) 

Here, h{p) = — ^ is the specific enthalpy, c{p) = y/P'{p) is sound speed, V = (01,82,83), and Co = 

-Qq + h{po) is Bernoulli's constant which is determined by the incoming uniform supersonic flow at minus 
inflnity with velocity (0, 0, go) and density po > 0. 

c^(p) 

By (1.1) and the implicit function theorem, in view of h'{p) = > for p > 0, the density p(x) can be 

expressed as 

p=h-^(co-^\V^\''^ =H{V^). (1.2) 

3 

Substituting (1.2) into the equation di{pui) = 0, which expresses the conservation of mass, yields 

i=l 

3 

J2diiH{V^)di^) =0. (1.3) 

i=l 

More intuitively, for any C^-solution (1.3) can be rewritten as a second-order quasilinear equation, 

3 

^((9i$)2-c2)92$ + 2 Yl di^dj^dfj^ = 0; (1.4) 

i=l l<i<j<3 

here c = c{p) = c(iJ(V$)). Note that (1.4) is strictly hyperbolic with respect to the a;3-direction in case 
M3 > c{p) holds. 

For the geometry of the conic surface, it is convenient to work in the cylindrical coordinates {z, r, 6), where 

a;i=rcos^, X2 = rsin^, X3 = z, (1-5) 
r = ^/x1+ X2, and < ^ < 27r. Under the change of coordinates (1.5), Eq. (1.4) becomes 

((5,$)2 - c2)52$ + {(dr^f - c')d^^ + 1 (^i^ - 0^)92$ + 2d,^(^dr^dl<^ + ^de^d'^e^ 

+ ^dM^d^e^ - -dr^ fi^^ + A = 0. (1.6) 



Let ^~{z, r, 9) and $+(2:, r, 6) denote the velocity potential for the flow ahead and past the resulting shock 

front r = x{z,d), respectively, where x(Oj ^) = 0- Then (1.6) splits into two equations. That is, ^^{z,r,0) 
satisfy the following equations in their corresponding domains, 

_ V-f^^^ + (c-)^Uo inn_ (1.7) 
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and 



V+f^^^ + (c+)^')=0 mn+; (1.8) 



here c± = c(iy(V$±)), Vt_ = {{z,r,e):r > x(z,6l),0 < 6 < 2-k,z > 0}, and 0+ = {{z,r,e):hQZ < r < 
x{z,e),Q< e<2Tx,z> 0}. 

On the conic surface r = boz, satisfies the boundary condition 

dr^+ - 6o9z$+ =0 on r = boz, (1.9) 

while on the conic shock F = {r = x(-2^)^)}) by Eq. (1.3) and the change of coordinates (1.5), the Rankine- 
Hugoniot condition becomes 

[H{V^)dr^]-[H{V^)d,^d,x=^[HiV^)de^]dex on T. (1.10) 
Moreover, the potential ^{z, r, 6) is continuous across the shock, i.e., 

^+{z,x{z,e),e) = ^-{z,x{.z,e),e) onr. (i.u) 

Furthermore, we impose initial conditions on '^~{z, r, 0), 

$-(0,r,^) = £$o(r,e), d,<^-{Q,r,e) = qo + e<^^{r,e), (1.12) 

where £ > is a small parameter, qo > c{po), and $(^(r, ^), $^(r, ^) € C^[{0,1) x [0, 27r]) for some fixed 
number Z > 0. 

The main result states: 

Theorem 1.1. For small ba > and a sufficiently large speed qo, there exists a small constant £o > 
depending on qo, po, bo, and j such that problem (1.7)-(1.8) together with (1.9)-(1.12) possesses a global C'°° 

supersonic shock solution {^^{z,r,6),x{z,9)) for any e < eq- Moreover, (Va;$+, ^ — ) approaches the 

corresponding quantities for the incoming uniform supersonic flow (0, 0, go) with density po past the sharp 
circular cone r = boz with rate (1 + z)~"^° for any positive number mo < 5. 

RemsLrk 1.1. TTie various smallness assumptions in Theorem 1.1 can be expressed by saying that 

< e < mini -^(6ogo)"^ ■4(?>o9o)"^ I < &o and 60 < 6*, 

where 6* is the critical value given in Remark 2.1 below. 

Remark 1.2. As in [5], [9-10], and [31], where suitable symmetry assumptions on the incoming supersonic 
flow with small perturbation or artificial boundary condition on the conic surface were imposed, we emphasize 
that also in the case treated in this paper there are no discontinuities for the weak solution $(x) = #"'"(x) in 
f2+ and <E>~(.t) in but the m,ain m,ulti-dim,ensional conic shock front. This m,eans that the supersonic conic 
shock is structurally stable in the whole space for a polytropic gas and arbitrary perturbations. This agrees with 
observations from physical experiment and numerical simulations. 

4 



Remark 1.3. The nonlinear hyperbolic equation (1.4) is actually a second-order quasilinear wave equation 
in two space dimensions if one regards X3 as time, as the flow is supersonic in x^-direction. By a direct 
verification, one sees that (1.4) does not fulfill the "null- condition" put forward in [6] and [19]. Therefore, 
in term,s of the results of [1], [12-13], [15-17], [25-26], and [28], if there was no shock for the solution to 
Eq. (1.4), then the classical solution to (1.4) would blow up in finite time. Thus, the result of Theorem 1.1 
asserts that the m,ulti- dimensional supersonic shock absorbs all possible compressions in the flow and prevents 
the flow from forming further shocks as well as other singularities. 

Remark 1.4. As energy estimates fail to hold in BY-spaces for multi- dimensional hyperbolic equations and 
systems as shown in [27], the method used in [21] {which was the Glimm scheme for a spherically symmetrically 
perturbed conic surface) cannot be applied to the genuinely multi- dimensional problem treated here. 

Remark 1.5. It was indicated in [7, pages 317-318 and 414] that if a supersonic steady flow approaches 
from minus infinity and hits a sharp cone in direction of its axis, then it follows from the Rankine-Hugoniot 
conditions and the physical entropy condition, by an application of the method of the apple curve (see Figure 3 
below) that there possibly occur a weak shock and a strong shock attached to the tip of the cone. These shocks 
are supersonic and transonic, respectively. It was frequently stated that the strong shock is unstable and that, 
therefore, only the weak shock is present in real situations. In [32-34], the global instability of an attached strong 
conic shock in the whole space was systematically studied {in this case, the corresponding subsonic potential 
equation is nonlinear elliptic and the steady Euler system becomes elliptic-hyperbolic) which especially showed 
that a global strong conic shock is actually unstable as long as the perturbation of the sharp circular cone 
satisfies som,e suitable assumptions. On the other hand, from the result in this paper, one infers the global 
stability of a supersonic conic shock. Consequently, in regard to the global stability or instability of weak and 
strong conic shocks, these results give a partial answer to the problem of Courant and Priedrichs [7] stated 
above. 











/>\ A' 






I stwr.k polar \ 








1 







Figure 3. Apple curve showing all possible end velocity of a conic shock. 

Remark 1.6. For large qo, the incoming flow is called hypersonic. The famous independence principle 
for large Mach numbers {that there exists a stable limit state for a hypersonic flow as the Mach number goes 
to infinity) is likewise illustrated by Theorem 1.1 for a hypersonic gas past a sharp cone. For more physical 
properties of hypersonic flow, see [2] , [8] , and [29] . 

Remark 1.7. In case the conic surface is also perturbed, where the multi- dimensional perturbation is small 
and of compact support {possibly including a compact perturbation near the vertex of the cone, as local existence 

5 



of such a conic shock was established in [4]) or decays svfficiently fast when z goes to infinity, a result analogous 
to Theorem 1.1 remains in force, where the proofs of this paper still work. 

Let us mention some work which is directly related to this paper. In [5] and [9], under the assumptions of 
an incoming uniform supersonic flow and that the angle of the spherically symmetrically curved conic body 
is sufBciently small (smaller than the critical angle which guarantees that the supersonic shock is attached), 
it was shown that a spherically symmetric supersonic conic shock exists globally past the conic body for a 
supersonic polytropic flow. Z. Xin and H. Yin established in [31] the global existence of a multi-dimensional 
supersonic conic shock for an incoming uniform supersonic flow past a generally curved sharp cone under an 
artificial boundary condition - the Dirichlet condition for the potential on the conic surface (physically, this 
kind of boundary condition means that the conic body is perforated or porous) . It should be emphasized here 
that the Dirichlet boundary condition for the potential in [31] played a crucial role in deriving a priori energy 
estimates and further obtaining the global existence. It means that the Poincare inequalities are available on 
the shock surface and in the interior of the downstream domain. By applying the Glimm scheme, in the case 
of a spherically symmetrically curved cone, W.-C. Lien and T.-P. Liu in [21] obtained the global existence 
of a weak solution and the long-distance asymptotic behavior under suitable restrictions on the large Mach 
number, the sharp vertex angle, and the shock strength. The main interest here is to establish the global 
existence of a genuinely multi-dimensional supersonic conic shock for a perturbed supersonic polytropic gas 
past an infinitely long cone with a sharp angle when the speed of the incoming flow is large. Especially, we 
remove the key assumption of spherically symmetry on the perturbed supersonic flow which was assumed in 
[5], [9], and [21] and which was essential in the proofs there. 

Let us also comment on the proof of Theorem 1.1. In order to prove Theorem 1.1, we intend to use continuous 
induction to establish a priori estimates of the solution and its derivatives. To achieve this objective, as in [5], 
[11], and [15], we need to derive global weighted energy estimates for the linearized problem (1.7)-(1.8) with 
(1.9)-(1.12). Based on such estimates, one then obtains the global existence, stability, and the asymptotic 
behavior of the shock solution to the perturbed nonlinear problem. The key ingredients in the analysis to 
obtain weighted energy estimates are an appropriate multiplier and a new Hardy-type inequality on the shock 
surface. Finding a suitable multiplier is a hard task for the following reasons: First, to obtain the global 
existence requires to establish global estimates, independent of z. of the potential function and its derivatives 
on the boundaries as well as in the interior of the downstream domain. This implies strict constraints on the 
multiplier and makes the computations delicate and involved. Secondly, as our background solution is self- 
similar in a downstream domain and strongly depends on the location of the boundary of the cone, the angle 
at the vertex of the cone, the Mach number of the incoming flow, and the equation of state of the gas under 
consideration, one needs to take some measures to simplify the coefficients of the nonlinear equation together 
with the corresponding nonlinear boundary conditions so that the procedure to find the multiplier becomes 
managable. Thirdly, for the multi-dimensional case, the Neumann-type boundary condition (1.9) fulfilled by 
introduces additional difficulties compared to [5] and [31], where [5] only treats the case of a spherically 
symmetric conic shock with Neumann-type boundary condition on the conic surface, while [31] treats the case 
of an artificial Dirichlet-type boundary condition for the potential on a multi-dimensionally perturbed conic 
surface. The latter plays a key role in the analysis of [31], as the corresponding Poincare inequality is available 
on the shock surface and the interior of the downstream domain, respectively, while this is not the case in the 
problem treated here. Thanks to some delicate analysis accompanied by a new Hardy-type inequality derived 
by making full use of the special structure of the shock boundarj^ conditions (i.e., the sizes as well as the signs 
of the coefficients in (3.6)-(3.7)), we finally overcome all these difficulties and obtain a uniform estimate of 
II Vx^"*" • From this, higher-order estimates of Vx^'*' can be established by using Klainerman vector field 

and commutator arguments together with a careful verification that suitably chosen higher-order derivatives 
of the solution satisfy the Neumann-type boundary condition on the conic surface. This eventually establishes 
Theorem 1.1. 

The paper is organized as follows: In §2, we derive some basic estimates of the background self-similar 
solution in case of an incoming hypersonic fiow, which are required to treat the linearization of the nonlinear 
problem and for the construction of the multiplier. In §3, we reformulate problem (1.7)-(1.12) by decomposing 
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its solution as a sum of the background solution and a small perturbation ip so that its linearization can be 
studied in a convenient way. In §4, we first establish a uniform weighted energy estimate for the corresponding 
linear problem, where also an appropriate multiplier is constructed. Based on such an energy estimate, we 
obtain a uniform weighted energy estimate of for the nonlinear problem through establishing a new 
Hardy-type inequality. In §5, by the estimates derived in §4, we continue to establish uniform higher-order 
weighted estimates of V^jt^. In §6, the proof of Theorem 1.1 is eventually completed by utilizing Sobolev's 
embedding theorem and continuous induction. Some lengthy computations are carried out in an appendix. 

In what follows, we will use the following conventions: 

. C stands for a generic positive constant which does not depend on any quantity except the adiabatic 
constant 7 (1 < 7 < 3). 

. C(-) represents a generic positive constant which depends on its argument (or arguments). 
. O(-) means that |0(-)| < C"! • | holds true. In particular, 0{e) abbreviates \0{e)\ < C{bo,qo)e. 
. dS stands for the surface measure in the corresponding surface integral. 



§2. Analysis of the self-similar background solution 



In this section, we will provide, with more details than in [5] , properties of the background solution when the 
Mach number of the incoming supersonic flow is large and the supersonic shock is attached. These properties 
will be applied again and again in the later analysis of §3-§5. 

Following the illustrations of [7, page 407] , the supersonic conic shock phenomenon for an incoming super- 
sonic flow past a sharp circular cone is described as follows: Suppose that there is a uniform supersonic flow 
(0, 0, (/o) with constant density po > which approaches from minus infinity. Let the flow hit the circular 
cone {(r, z): r < h^z, 2; > 0} in the direction of the z-axis. Then there exists a critical value h* > 0, which is 
determined by the parameters of the incoming flow, such that there occurs a supersonic conic shock r = .sqz, 
.So > fco; attached to the tip of the cone whenever 60 < b* holds true. Moreover, the solution to (1.3) with 
(1.1) past the shock surface is self-similar, that is, in cylindrical coordinates {z,r,9), the density and velocity 
between the shock front and the conic surface are of the form p = p{s), ui = Ur{s) — , ui = Ur(s)— ■) and 



"3 = Uz{s), where s 
differential system: 



In this case, Eq. (1.3) with (1.1) can be reduced to the following nonlinear ordinary 



p'{s) - 
<(s) 



pUr{sUz — Ur) 



S((l + S2)c2(p)-(SM,-M,)2)' 

C^{p)Ur 

s{{l + S^)c\p)-{sU,-UrYy 

C^{p)Ur 

(1 + S2)c2(p)-(SU,-U,)2' 



for bQ < s < Sq. 



(2.1) 



As explained in [5] or [7], for the denominator it holds (1 -|- s^)c^{p) — [suz — u^Y > for 60 < s < sq, which 
implies that system (2.1) makes sense. 

On the shock front r = sqz, it follows from the Rankine-Hugoniot conditions and Lax's geometric entropy 
conditions on the 2-shock that 

I [mz]+soK]=0 



and 

( Ai(so) < -So < A2(so), 



where 



Additionally, the flow satisfies the condition 

Uribo) = bou^{bo) (2.5) 

on the fixed boundary s = bo. 

As indicated in [7, pages 411-414] or [18], the boundary value problem (2.1)-(2.4) can be solved by the 

shooting method as well as by the method of the apple curve. More specifically, for any given ba > 0, which is 
smaller than the critical value 6*. one can determine the solution to (2.1)-(2.4) by finding the integral curve of 

— — = from the intersection point of the apple curve with the ray u,- = bgUz to some point at the shock 

dUz Ur 

polar (see Figure 3 above or [7, Fig. 8 on page 414]). In this paper, such a supersonic solution past the shock 
is called the background solution. 

For large qo, some detailed properties of the background solution are as follows: 



/I /7 + 7 

Lemma 2.1. For qo large enough, 1 < 7 < 3, and < 60 < ^* = \ 7;{\ ; — 1); one has, for 



2 



bo < s < So, 

(i) so = bo{ 1 + O((6o9o)"^) 



ii) < suzis) - Uris) < 0((6oQo)^). 



iv) «.(«) = Yf^(l + 0((Mo)- — ) 



''^'^ = ( 2^7(1 + 6g) ) (^0*)^ + Oiibogo)-'') + Oiiboqo)-^) 

vi) q\8) - cHpis)) = ql ((1 + O{{boqo)-^) + O((6o9o)-^)) , 

here and below (f{s) = u^ls) + u1{s). 



vii) ulis) - cP{p{s)) = L_!l!^^g±Ml^2 ^ 0((6ogo)-2) + O{{boqo)-^) ) > 0. 

7-1 



(1 + 6, 



viii) (1 + s'y{p{s)) - {suz{s) - Ur{s)y = ^{boqoY ( 1 + Oiiboqo)-') + 0{(boqo)- — ) ) > 0. 



Remark 2.1. ^* ~ \j 2^y ^ ^ — 1) is a root of the quartic equation 1 — ^L_L6q(1+6q) = 0. This guarantees 

that the flow across the shock is still supersonic in direction of z for bo < 6* and large qo ■ It can be derived 
from the expression in Lemma 2.1 (vii). 

p+ 

Proof. Set p+ = lim p{s), m+ = lim Ur{s), m+ = lim Uz{s), and a = — . It follows the Rankine- 

s— >so— s— >so— s—^sq— Po 

Hugoniot conditions (2.2) that 
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'^'^) (2.6) 
+ soqo L. 1 



Substituting (2.6) into Bernoulli's law (1.1) yields 



ip^r^' - pr\p^r + tttt^ pt - (p^r = o. (2.7) 



7-1- 7-1"" ' 2(1 + 4) 



7-1 V ' V. . y . 2(l + s§) 

Set 



Then (2.7) implies that F{po) = F{p'^) = 0. It follows from a direct computation that, for x G (0, po), 



-47(7+1) . 2A7 1 52^2 / A7(7 + l) 2^7 



< J^^pr - ^pV - ^) = 4c\p.) ''"'^ 



V 7-1 7-1 1 + So/ V 1 + s 

Due to Lax's geometric entropy conditions (2.3), F'{x) < for a; G (0, po). Together with F{po) = F{p'^) = 0, 
this yields /?"*" > po- 

In this case, (2.7) is equivalent to 

, a^-^-l _ {^-l)slql 
" a^-1 - 2^7(1 + .^)pr' ^'-'^ 

where a > 1. Since the left hand side of (2.8) is bounded if a > 1 is bounded, for go is large, a is also large. 
From this fact, one has 

" = ^ (^I^^TR)) ~^ + OiiM.)-') + 0((Mo)-^) 

Substituting this into (2.6) yields 

l + 0((6o(Zo)"^) 



•sogo 



(2.9) 



Moreover, from u'^{s) < and u'^{s) > for s G [bo, sq], one has 

Ur < Ur{s) < Ur{bo) = boUz{bo) < &OMz(s) < bouf . (2-10) 

Combining (2.9) with (2.10) yields (i) and further (iii)-(iv). 
(ii) comes from (2.9) and the fact that 

= boUz{bo) — Ur{bo) < suz{s) — Ur{s) < squ'^ — for bQ < s < So; 

here {suz{s) — Ur{s))' > Q ior bo < s < sq has been applied. 

Furthermore, Bernoulli's law (1.1) and (i)-(iv) show by a direct computation that (v)-(viii) hold. Therefore, 
the proof of Lemma 2.1 is complete. □ 

Lemma 2.2. Under the assumptions of Lemma 2.1, one has, for bo < s < sq, 

r^w^ VY^jl + bl)bo{VY^bl - ^2 - - l)b'',) , ^^^^ ^-2^^nrrA ^-^^^^n 
(i)Ai(s)-s= 2- (7- 1)62(1 + 62) (1 + O((6ogo) ) + O((6ogo) 7-1)) <0. 
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(ii)A2(s)-s = ^_ -^.^^2(;^^^2-) '-[1 + O{{boqo) )+0((6o9o) ^"Oj > 0. 



(iii) = -(^^^2)2(1 + Oiiboqo)-') + Oiiboqo)-^- 

(iv) = 7T-^fl + 0((Mo)-')+0((6ogo)-^) 



+ 

(v) |p'(s)|<0(i). 

Proof, (i)-(v) can be directly derived from (2.1), (2.4), and Lemma 2.1. For the reader's convenience, we 
provide the detailed computation for (i) as an example. 

It follows from (2.4), (iii)-(v) in Lemma 2.1, and a direct computation that 

Uris)u,is) - C{p{s))y/U^{s) + ul(s) - c2(p(s)) 



bpgo go 

1+6: 



Ulis) - c2(p(,s)) 

7/^ 



'l + 0((6o«o) ) + C((6o9o) -'"O 



2bo - (1 + bg)Vy^bo V2 - (7 - l)b] 
2-(7-l)6§(l + 6g) 



O((6o(zo)'') + O((6o9o)"~) 



Due to So = &o + 0{{boqo) ) j by (i) in Lemma 2.1, Ai(s) — s satisfies (i). □ 

Remark 2.2. 5mce the denominator of system (2.1) is positive in the interval [bo, sq] (see Lemma 2.1 (viii)), 

one can extend the background solution (p{s),Uz{s),Ur{s)) to (2.1)-(2.3) and (2.5) to the interval [60, sq + tq] 

_ 4 

for some small positive constant tq satisfying < To < Qo (■'^0 — bo). In the following sections, we will denote 

this extension of the background solution to {(z, r):z > 0, boz <r< (so + to)z} by (p(s), 'U^(s), 'Ur(s)), where 
r 

s = - . The corresponding extension of the potential will be denoted by $(s). 

§3. Reformulation of the nonlinear problem 

In this section, we reformulate problem (1.8)-(1.11) by decomposing its solution as a sum of the background 

solution and a small perturbation. Moreover, based on the analysis of the background solution in Lemmas 2.1 
and 2.2, we establish estimates of the coefficients which appear in the reformulated problem when qo is large. 

We now provide a global existence result for the solution to Eq. (1.7) with initial data (1.12) ahead of the 
shock. 

Lemma 3.1. Eq. (1.7) with (1-12) possesses a C°° -solution <^~{z,r,6) in Jl_ — {(z, r, 0):r > x(2;, ^), < 
< 27r, z > Q). Moreover, ^~iz, r, 6*) — qoz G C^(f2_) and, for any fc G N, there exists a positive constant Ck 
independent of e such that 

||$-(z,r,6')-go^||c'=(f2_) < Cue. 



Proof. The quasilinear equation (1.7) is strictly hyperbolic with respect to the direction of z, for the 
supersonic flow condition Ug > c{p~). The initial condition (1.12) means a small perturbation with compact 
support away from the origin. Thus, for large qo, by finite propagation speed and standard Picard iteration 
(or see [16]), one easily derives the validity of Lemma 3.1. □ 

Note that there exists a constant Tq > such that {z,r,9) = qoZ ioT z > Tq. Without loss of generality, 
below we will assume that To = 1. 

Next, we reformulate the nonlinear problem (1.8)-(1.11). For notational convenience, we will neglect all the 
superscripts "+" in (1.8)-(1.11) from now on. 
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Let $ be the solution to (1.8)-(1.11) and <p be the perturbation of the background solution, that is, = 
here $ is given in Remark 2.2. Then, by a direct computation, (1.7) is reduced to: 



where 



and 



It* i 1 1 'T* 

+ -/5(-,V:,(/))9fe<^+ -/6(-, V^v9)a^eiy3+ -/7(-, Va;(p) in 

r Z T Z T z 



P2{S) 

Psis) 

P5{S) 



Ul{s))-C^{p{s)y 

ui{s) - c^{p{s)y 

1 / 7 + 1 



Ul{s) - C2(p(s)) 



7 ~ 1 

2 s Uz{s)u'^{s) + —^su;,{s)u'^{s) 



7-1 

+ SUr{s)u'^{s) ^ ^ — Uz{s)Ur{s) 



1 



ul{s) - c?{p{s)) 



-'^^-^s'^Ur{s)u'^{s) + '^^-^SUr{s)u'^{s) 



+ suz{s)u'z{s) + ^y-!-u^(s) - ]^c^{p{s)) 



/l(s, Va;iy?) 

feis, Vx'p) 



ul{s) 



-c^{p{s)){ 



7-1 



{2uz{s) + dzifi)dzip 



1 [ ... .«,.^2 7-1 



-2Ur{s)drip — {dr<pY 



{2uzis) + dz<f)dzip 



{do9? ,7-1 



1 



ul{s) - c^{p{s)) \ r 



+ 



{2uz{s) + dz(p)dz(p 



ulis) - c^pis)) 
2 

U2(s) - C2(p(s)) 



{I 



Uz{s)de(p+ -dz(pd0(p 
r 
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(3.1) 



(3.2) 



where s = -. 

z 

On r = boz, one has 

^riyj = bodzip. (3.3) 
On the free boundary r = x(z, by the continuity condition (1.11), (1-10) can be written as 

Ji-(V$)((a,$)2 + (5,$)2 - qod,'^) - poqod.^ + poql = - {V<^>){de^f on r = x(^,^). (3.4) 

Introducing the notation 

z 

Then (3.4) can be rewritten as 

Bi9^<yj + B29^9' + S3C = K(C,Va.</j) onr = x{z,e), (3.5) 

where 

^1 = 57^7^fwr(so) + Wz(so)(m2(so) - Qo)V'-('^o) + 2p(so)Mr(so), 



c^(p(«o)) 

-B2 = — ("r(go) +M2(so)(mz(so) - 90) )m2(so) + 2p(so) (Mz (sq) - 9o) + (K^o) " ^0)50, 

B3 = p{sq) ^2wr(soX(so) + 2(u2(so) - Qo)uz'{sq) + goWz'(so)^ + p'(so) ^Wr(so) + Wz(so)(w2(so) - Qo) 
- PoqoUz'iso), 

and the generic function k(^, Vx^y?) is used to denote a quantity dominated by C(6o, ?o)|(Ci '^x'fi)\'^- 
By Lemma 3.3 below, one has Bi ^ in (3.5) for large qq. Thus, Eq. (3.5) can be rewritten as 

Bo(p + IJ.2^ = KiCVx^p) on r = x{z,0); (3.6) 

here Bo(p = df^p + p-id^ip with /ii = — and /i2 = — . 

Bi Bi 

Furthermore, (1.11) and the properties of $~ for z > 1 and 6 G [0, 27r] listed in Lemma 3.1 imply 

ip{z,x{z, o),e) = ^z, xiz, e),e)- s^z, e) - (4(z, x(^, e), e) - ^z, s^z, e)) 

= -{^z,x{z,O),0)-Hz,soz,e)) 
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On the other hand, it follows from Lemma 3.1 that — qqz G Cq°{Q,-) in (1.7) and that, near the vertex of 
the cone r = boz, the solution $+(x) is actually the background solution. Thus, in order to prove Theorem 1.1, 
by the local existence and stability result for the mult i- dimensional shock solution to the potential flow equation 
in [24], we only need to solve problem (3.1) in the domain {(z, r, 9): z> 1, boz < r < x{Z; 0), < 9 < 27r} with 
the boundary conditions (3.3), (3.6)-(3.7), and small initial data (p{z, r, 9)1^=1, dz^p{z, r, 9)\z=i, and ^{z, 6')lz=i. 
Here, smallness of initial data means that 

J2 I VL9^(l, r,9)\ + J2 I VU(1' ^)l < (3-8) 

l<ko l<ko 

where fco e N, fco > 7. 

For later use, wc now list specific estimates of the coefficients in (3.1), (3.2), and (3.6). Since these estimates 
result from a direct, but tedious computation that makes use of Lemmas 2.1 and 2.2, we postpone the proof 
to the appendix. 

With respect to the coefficients of C<p in (3.2), one has: 

Lemma 3.2. For qo large enough, 1 < 7 < 3, < 6o < o,nd bo < s < sq, one has 

h) j , , ^^^^ „ ^-2^ 

1_ 2^52(1 + 62) 



= . .-1,,.. , ,,J i + o((Mo)-^)+o((Mo)"^) , 



''^'^ ^ (1 + ^((^0.0)-) + o((Mo)-- )) , 

^3(^) ^ ^^^1^ (1 + 0((Mo)-) + 0((Mo)-^ 
P4(s) = 0((Mo)"') + 0((6osor*). 

Moreover, Bi, i = l,2, 3), in (3.5) and jij, j = 1, 2, in (3.6) admit the following estimates: 
Lemma 3.3. For large qo, one has 

1 

= {^mrm) ^ ^'"^"^^ + """'"'"^"^ + omqo)-^)) > 0, 

^2 = , ]~!'L f ./V.^D (&o9o)^ f 1 + 0((foo9o)-') + ©((Mo)-^)") , 



60(1 + 6^) V2^7(l + &o) 

^3 = -rrr^-To^ f oiV.^L ) (&o9o)^ f i + i^O((Mo)-') + ^0((Mo)-^) ) , 



60(1 + 6^)2 V 2^7(1 + &o)/ V bo ' ' 60 

Ml = (1 + 0((&o9o)-') + 0((6o«o)-^)) > 0, 

= -2(TT6l) + l^iM-^) + lo((Mo)-^)) < 0. 
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Remark 3.1. We emphasize that there is a large factor qo in the coefficient /LI2 of the boundary condition 
(3.6) which shows that one has to be very careful in treating the shock boundary condition later on. 



§4. A first-order weighted energy estimate 

In this section, we establish a weighted energy estimate of Vx'P for the hnear part of (3.1), together with 
(3.3) and (3.6)-(3.8), which will play a fundamental role in our subsequent analysis. 

Set Dt = {{z,r,e):l < z < T, boz < r < x(^,^),0 < 6 < 2tt} for any T > 1. Ft = {(z,r,6'):r = 
X{z, d),l<z<T,0<d<2'K} and Bt = {{z, r,9):r = boz, 1 < ^ < T, < 6* < 27r} are the lateral boundaries 
of Dt. 

Theorem 4.1. Let (p G C°°{Dt) satisfy the boundary condition (3.3). Let \^{z,0)\ + \zdz(,{z,0)\ + 
\d0^{z,0)\ + |Vx<^| < Me hold true for small s, z £ [1,^], and some positive constant M. Then there ex- 
ists a multiplier M.(p = A{z,r)dz(p + B{z,r)dj4> with smooth coefficients such that, for any fixed constant 
H < -1, 



CiT^+i/"/" {Vxipf{T,r,e)drde + C2 Iff z^iV^^p)^ dzdrdO 

J JboT<r<x(T,e) JJJdt 

+ CsJJ^ z^+\dz^rdS + C,JJ^ z'^+'^^^dS 
< [[[ C^- M^dzdrdO + C5 ff ({dzv^ f + {dr'fi f + 

JJJDt J Jbo<r<x(l,e)\ 



{l,r,e)drde 



+ Ce z^^+^Boipf dS, (4.1) 
J Jtt 

where Ci, 1 < i < 6, are some positive constants depending on bo and qq. In particular, 



(7 - 1)6(1(1 + b^) 

2(1-2^62(1 + 62)) 



(4.2) 



7. + Oiiboqo)-') + 0((Mo)-^). 



Remark 4.1. The values of constants C3 and Cg will play an im,portant role in the energy estim,ates for the 
nonlinear problem (3.3) with (3.6)-(3.8). The most troublesome term J jp z^'^^ [Boip)'^ dS in the right-hand side 
of (4.1) will be shown to be absorbed by positive integrals in the left-hand side of (4.1). The reason for which 
the term JJ^ z^^~^^ {Boip)"^ dS is hard is the following one: In view of the Neumann boundary condition (3.3), 
which is different from the artificial Dirichlet boundary condition used in [31], the usual Poincare inequality 
is not available for the solution (p on the shock surface which means that the L'^{Tt)- estimates of'Vx'P cannot 
be obtained directly. (Note that the boundary condition (3.6) contains the function ^, and that ^ is roughly 

equivalent to — in view of (3.7), so that an estimate of ip on the shock surface has to be established.) 
z 

Proof. We will determine the coefficients A{z, r) = z'^ra{j) and B{z, r) = z'^'^^b{-) for z >1 later on. Set 

M.if = A{z, r)dzip + B{z, r)dr^p. Then it follows from an integration by parts that 

JJJ Cip-Mipdzdrde = jjj (^-i^{dz^pf + K2dz<pdr<p + K^idr^pf + K4, ^^'^^J ^ dzdrdO 
+ T''+^ ff Ni{T,r,e)dS- ff Ni{l,r,e)dS 

■' ■'b„T<r<x{T,e) J Jbo<r<x{1.9) 



+ 



JJ^ z^+i (^N2 - dzxNi - dexN^ dS + jj^ z^'+^ (^boN, - N^^ dS, (4.3) 
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where 



Ki =[^s^- sPiis) )a'{s) + + -f s - Pi(s) - sPi{s) + 2Pi{s) ]a{s), 



K2 = -sP2is)a'{s) + sb'is) + (^-P2(s) + 2P5(s) - sPi{s)ja{s) + (^-(m + 1) + f P4(s)j6(s), 

Ks = -ls^P2{s)a'is) + (^sP.is) - iP2(s))6'(s) + (^f sPsC^) - is2p^(s))a(s) (4.4) 



and 



-ill + mis) + sPiis) - \p^{s) + ip,is)jbis), 

K, = ls^Psis)a'{s) - iPsisWis) + (-f sPgC.) + is2p^(.))a(s) + (^^Ps{s) - iP^(«))6(«) 

Ni = isa(s)(9,<^)2 + bis)d,^dr^ + f 6(s)Pi(s) - hais)P2is)) {dr^f + sais)P3is) 



2r2 ' 

iV2 = [sa{s)PM - \h{s)^ {d^^f + sa{s)P2{s)d,^dr^ + i6(s)P2(s)(a.^)2 + 6(s)P3(s)^^, ^^-^^ 

-^Va = -^sa{s)P3{s)dzipde'f> - ^b{s)Pz{s)dr'pdeip. 

Our goal is to choose positive functions a{s) and 6(s) in such a way that all integrals over Dt, Bt, and 
z = r in the right-hand side of (4.3) are either positive or zero and that the integral over Ft provides sufficient 
control on (p. Starting from here, we will derive sufficient conditions on the choice of o(s) and b{s) in the process 
of analyzing each integral from which o(s) and b{s) can then be determined. This process is subdivided into 
the following five steps. 

Step 1. Handling of J J z'*"'"^ ^&o-^i — -^2^ dS. By the boundary condition (3.3) and Ur{bo) = boUz{bo), 
one obtains that 

, .2s , bl{l + bl)ulibo) 

'^'°>^2iul{bo)-c\p{bo))). 



bom -N2={-Jl + bl)+ ) (bibo) - bla{bo)){d.^)' 



-(6(6o)-6ga(6o))P3(M^^ 

on Bt. Since the terms | and P3 (60) are positive according to (iv) and (vii) in Lemma 2.1 

2(u^(6o) - c^(p(&o))) 

and Lemma 3.2, in order to control the integral J J z^~^^ (^oNi — dS, one should choose 

b{bo) = blaibo). (4.6) 

This then implies 

J I ^"+1 {boN^ - dS = 0. (4.7) 
Step 2. Positivity of / / N^dS. To ensure the positivity of the terms in A^i, which are 

J JboT<.r<,x{T,e) 

quadratic in {dz(p,dr<fi, — ^), a(s) and 6(s) should fulfill 

r ais) > 0, 

1 b^is) - 2sPi(s)a(s)6(s) + s^P2is)a^is) < 
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which is equivalent to 



In this case, one arrives at 



a{s) > 0, Ai(s) < < A2(s). (4.8) 

sa{sj 



[f N,dS> [[ (xrnUs){{dzV'f + {dr^f) + ^saP3{s){de^)AdS, (4.9) 

J JboT<r<x(T,e) J J boT<r<x(T,e) \ ^'^ / 

where 



Xminis) = \ Qsa(s) + 6(s)Pi(s) - \sa{s)P-,{s) - ^J i^sais) - 6(s)Pi(s) + \sa{s)P2is)f + b{sy^ . 

Moreover, by Lemmas 2.1, 2.2, and 3.2, the assumption on ^{z, 9) in Theorem 4.1, and the choices of a(s), 6(s) 
to be made later on, one actually obtains that 

Amin(s) >C7(6o,9o)+0(£). (4.10) 

Step 3. Positivity of the integral over Dt- Under the constraints (4.6) and (4.8), we will choose o(s) 
and b{s) in such a way that 



This will be true if the coefficients Ki, 1 < i < 4, satisfy 



Ki >0, 

- ^iKs < 0, (4.11) 
Ki > 0. 



According to (4.6), one can choose 

o(s) = 1, 



b{s) = s^l + -{s-bo) 



where the constant b will be determined in a moment. 

It follows from Lemmas 2.1, 2.2, and 3.2 and a direct computation that 



(4.12) 



Ki = (^-6o + 0{bt)j & - (f - l)&o + 0{bi) + 0{{bm)-') + O((6o9o)-^), 
K2 = (bl + 0(63)^6+ (2 - ^l)bl + 0{bl) + O((6oao)-') + O((6o<zo)"^), 



^3 = ( ^^hl + 0{bt)\ 6 + f 1 - M2±ll) bl + 0{bt) + O((6o9o)-^) + O((6o?o)-^), (4-13) 



Ki = - {^-^bl + 0(6^)^ 6 - ^(7 - l)bl + 0(6^) + O((6ogo)-') + ©((Mo)"^), 

- 4JCiif3 = -^(m - &)(m - 2 - 6)6^ + 0(6^) + O((6o9o)-') + ©((fcogo)"^). 
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Therefore, one obtains from (4.11) that, for qg large enough, 

6 + 2-/i>0, {fi-b){ii-2-b) > 0, 6 + /x<0. 

If one sets 

6=1, (4.14) 

then 

At < -1. (4.15) 

In this case, one arrives at 

> Cih) J ({dz^f + [dr^f + dzdrdO. (4.16) 

Step 4. Estimate of J J ^ dzX^i — dgxN^ dS. By the assumptions on ^{z, 6) in Theorem 

4.1, it follows from the expressions for N^, N2, N3 and a direct computation that 

N2 - dzxNi - dexNs = Mdz'fi? + Pidz<fdr<p + p2{dr<fi? + Pz^{de'ff, (4.17) 

where 

Po = soPi(so)a(so) - ^So«(«o) - ^6(so) + 0{e), 
/3i = soP2(so)a(so) - sob{so) + 0{e), 

132 = ^P2{so)b{so) - soPi{so)b{so) + isgP2(so)a(so) + 0(e), 
^3 = \P3{soMso) - sla{so)) + 0{e). 
Due to drip = Bo(p — fJ^idz^, from (4.17) one obtains that 

N2 - dzxNi - dexNz = {Po - Mi A + p2f4){dz^f + (A " 2f,ip2)dz^Bo^ + p2{Bo^f + Ps^ideipf. 
By Lemmas 2.1 and 3.2 and the expressions for a(s), b{s) in (4.12) and (4.14), one arrives at 

Po - f^iPi + P2fii = ^.^ ^[^.^ , i^blil + blf + O{{boqo)-^) + O((6o«o)-^) + 0{e] 
4(1- V^o(l + &o)) V 2 

Pi - 2/ii^2 = O((6o9o)"') + O((6o?o)"^) + 0(e), 

1 f 7-1,4/1 , u2\ , ^//j, „ ^-2^ 



^2 = i_ 7^52(1 ^^,2) ( -^^0(1 + bl) + O((6o9o)-') + Oiiboqo)-—^) + 0(£) ) , 

p3>C{bo,qo) + Oie). 
Consequently, 



^'^+1 (^N2 - dzxNi - dexNs^ dS 



>{Qi{bo,qo) + 0{e)) // z'^+\d,^fdS-{Q2{bo,qo) + 0{e)) // z''+\Bo^f dS 

J J ■jp J JT^'j' 

+ (0(60, 90) + 0{e)) j j z^+'^^dS, (4.18) 
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where 



Qi(6o,?o) 
Q2ibo,qo) 



(^_l)b§(l + b2)3 

8(1-2^62(1 + 62)) 

(7- 1)6^(1 + 6g) 
2(1-2^62(1 + ^,2)) 



+ O((&o?o)"') + O((&o9o)"^) > 0, 
+ O((&ogo)"') + O((6ogo)"^) > 0. 



ion 



Step 5. Estimate of / / Ni{l,r,e)dS. iProm the expression for Ni and the initial conditi 

J Jbo<r<xii,0) 

(3.8), one easily obtains 

// K5{l,r,e)drde <C{bo,qo)e^. (4.19) 

J Jbo<r<x(lfi) 

Finally, substituting the estimates (4.7), (4.9)-(4.10), (4.16), and (4.18)-(4.19) into (4.3), (4.1) and (4.2) 
are obtained in terms of the smallness of e given in Remark 1.1. Therefore, Theorem 4.1 is proved. □ 

Based on Theorem 4.1, we will derive a first-order uniform energy estimate of \7xV' fo^" the linear part of 
(3.1) together with (3.3) and (3.6)-(3.8). To this end, we require a Hardy-type inequality on //p^ z'^~^\(p\'^dS, 
which is motivated by [14, Theorem 330] and derived utilizing the special structures of (3.6)-(3.7). 

Lemma 4.2. (Hardy-type inequality) Under the assumptions of Theorem 4.1, for /i < —1, one has 

h2\2 



JJ^ z^->pd5<C(6o,go)£^+ (l + 0((Mo)-^) + 0((&ogo)"^)) JJ^ z^+\d,ipfdS 



(4.20) 



Proof. First we assert that 



jj^ z'^-^\^\^dS<{^+0{e)) jj^ z'^+\l-iiAxf{d.^fdS 
+ C{bo,qo)e(^JJ^ z'^+\BoipfdS + jj^ 



,M+i^d5)+C(6o,?o)e^ 



Indeed, note that 



// z>'-^ip'^dS= [ d0 [ z''-^ifi'^{z,xiz,0),e)dz. 
J Jtt Jo Ji 



(4.21) 
(4.22) 



Set m{9) = z^ ^^'^{z, x{z, 0), 0)dz. Then, by an integration by parts 



mie) = -z'^ip\z,x{z,e),0) 

A* 



< 



,e)dz 



- - [ zi^ifiz, x{z, e), e){d,<fi + d,xdr<fi)iz, x{z, e), 
1 M Ji 

p^2(i^ ^(1^ ^)^ e)-^£ z^^iz, xiz, 0), 9) {d.xBo^ + (1 - ^ll^,x)^^^) {z, xiz, e),e)dz. 
By (3.7), one has 

^{z,xiz,e),e) = ai{z,x{z,e),e)^, 

where ^ 

ai{z, xiz, e),e) = -(^j^ uriso + t^z, e))dt^ z<o. 



(4.23) 
(4.24) 
(4.25) 



Due to the assumptions in Theorem 4.1 and Lemma 2.1, dzX = ^o(l + O((6o?o) + 0(e)), and thus 

because of the smallness of e, 

1 - m^.X = ^^(1 + O((6o9o)"^) > 0. (4.26) 
It follows from (3.8), (3.6), /i2 < by Lemma 3.3, (4.24)-(4.25), and a direct computation that 

m(e) < C{bo, qo)e^ - - / z'^'fiiz, xiz, e), e){l - n^d,x)dMz, X{z, 0), 9) dz 
M Ji 

+ ^ [ z'^dzX'fiz, x{z, 0),O)^dz - - / z^dzXV^iz, xiz, 0),O)k{^, V(p) dz 
= C{bo, qo)e^ - - I z^'fiz, xiz, e), e)il - fiid,x)dMz, xiz, 0), 6) dz 

+ ^/ z^'d,xaiiz,xiz,e),e)edz-- j z''d,x^iz,xiz,e),e)Kii,v^)dz 

< Ciho, qo)e' + z"- V'(^, xiz, 0), 0)dz+^ £ z^+\l - fi,d,xfid,<fifiz, xiz, 6), 6) dz 
+ C(6o, qo)e z>^-^ip\z, xiz, Q),e) dz + z^+^e dz + £ ^^+i(|V.,,0p 



^-^)iz,xiz,0),e)dz 



Together with (4.22) and (4.26), this yields 



<Cibo,qo)s' + i^ + Oie)) JJ^ z^+\l-,i^d,x)id.^?dS 
+ Cibo,qo)e( [[ z^+\Bo^rdS+ [[ z>^+^i^dS^ 

< C(bo,go)£' + ^f^^ (1 + 0((Mo)"') + 0((&ogo)"^)) // z''+\d,^fdS 

+ Cibo,qo)e^ll^ z^+\Bo^rdS + z^+^i^dSy 

Hence, Lemma 4.2 is proved. □ 

Theorem 4.3. Under the assumptions of Theorem 4.1, for fi < —1, one has 



CoT''+^ [ [ I Vx<p(T, r, e)\'^drde + Co [[[ z^'lV^if]^ dzdrdO 

J JboT<r<x{T,e) JJJdt 

+ Co[[ z^+^\V ^ip\^ dS < f f f £<fi-M<pdzdrde + Cibo,qoy, (4.27) 

J JVt J J J Dt 



' J J J D-j' 

where Co = Co(6o) 9o) > Q is a generic constant. 



Proof. To obtain (4.27), we are required to give a delicate estimate of the term Cq z^~^^iBoip)^ dS in 

J Jtt 

the right-hand side of (4.1) so that it can be absorbed by the positive terms in the left-hand side of (4.1). 
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We now treat the term / / ^"+^(^01^?)^ dS. 

J Jtt 

^Prom (3.6) and the definition of k{^, Vx'p), one has 

'(60,90) // Z^+^K^(X,^x^)dS. 

(4.28) 



<HU1 + O{{boqo)-^)) ff z^'-\z^fdS + Cibo,qo)[[ z''+^K\^,Vxip)dS. 



Note that 



I f z^'+^K\^,Va.^)dS < C{bo,qo)e^ [ [ z''+\e + |V.(^|^) dS 

< Cibo, qo)e' f f z^+' (e + id.^r + + (Bo^r) dS. 

JJft \ ^ / (4.29) 

Furthermore, as e is small and by the boundary condition (3.7) together with Lemma 2.1 (iii), one obtains 
from (4.28)-(4.29) that 

z'^+\B,<pf dS < ^'(^^^+|'p' (^l+0((bogo)-^) + 0((Mo)-^)) jj^ z'^-^^-'dS 
+ C{bo,qo)e'll^ ^.+i(^(5,<p)2+(^^dS 
^ 4 (1 + ^O((^o.o)-) + i^O((6o,o)-^)) / dS 

+ Cibo, qo)e' I z>^+' [{d,^f + dS. (4.30) 

Therefore, combining (4.30) with (4.20) in Lemma 4.2 yields 
// z'^+^Bo^fdS 

J JTt 

<C(6o,go)e' + ^^^5^(l + ^0((6ogo)-') + ^0((6o9o)-^)) jj^ z^+\d,^f dS 

+ Cibo,qo)ell^ zf^+'^-^dS. (4.31) 
Substituting (4.31) into (4.1), one obtains 



CoT^+i [ [ \Vx<ff drdO + Co fjf ^"1 Vx<^|^ dzdrdO 

J JboT<r<x(T.e) JJJdt 

+ (Qo(&o,«o) + ^0((Mo)-') + ^O((6o«o)-^)) JJ^ z^+\d,ipfdS 

+ C0JI z^^^^^dS 

< [[[ C(p-M(pdx + C{bo,qo)e'^, (4.32) 
J J Jdt 



IDt 
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where Qo(bo,<lo) = —, nT~^ ° „ . (l 5-) and Co = Co(6o5 9o) > is a generic constant. Moreover, 

8(1- V^o(l + ^o)) 

Qo{bo,qo) > 

because of > 1. Together with (4.32), this finishes the proof of Theorem 4.3. □ 

§5. Higher-order weighted energy estimates 

In this section, we derive a higher-order energy estimate so that one can estabhsh the decay properties of 
Vx'P and ^ for large z. 
As in [20], we denote by 

S= {81,82}, where 8i = zd, + rdr, 82 = de, (5.1) 

certain Klainerman vector fields and by 

•S'r = {•S'lr, S'2r}, where 8ir = zd^ + zdzx{z,0)dr, 82T = de -\- deX^r, (5.2) 

vector fields which are tangential to F. We will use these vector fields, which are tangential to the cone surface 
and are nearly tangential to the shock surface, respectively, to act upon Eq. (3.1) and the boundary conditions 
(3.3) and (3.6). This allows us to raise the order of the energy estimates by a commutator argument. Let us 
point out that there is a difference to the usual commutator argument inasmuch as the vector fields S are only 
nearly tangential to the shock front which causes certain error terms to occur in the estimates which in turn 
is due to the perturbation of the shock surface with r = sqz. Furthermore, we cannot adapt the analysis of 
[31] as we have to deal with a Neumann-type boundary condition on the fixed boundary, while [31] treats an 
artificial Dirichlet-type boundary condition so that there a Poincare-type inequality (see also [11, Lemma 1]) 
as one of the key elements of the analysis of [31] is available. However, by making use of the delicate energy 
estimate established in §4, we will be able to drive the desired estimates. 

To prove Theorem 1.1, we require the following elementary estimate: 
Lemma 5.1. Let <p be a C'^°{Dt) solution to (3.1), where ko G N, and 

o<;i+i2<[^]+i o<(<[^]+i 
where M > is some constant and e is sufficiently small. Then 

c E IVxSVl< E 4lvL+Vl<c(6o,go) E l^-^'^l i^^^- (^-3) 

o<;<fco-i o<;<feo-i o<;<fco-i 



Proof. The first inequality is immediate from the definition of the Klainerman vector fields 8 in (5.1). 
We show the second inequality in (5.3). The case ko = I can be verified directly. Cases when ko > 3 can 
then be obtained by an inductive argument. Thus, we only need to deal with the case ko = 2. 
Because of 

dr = ^ (^81 - zd;,^ , 

dl = ^(^d,Si-d,-zd^y 

d^ = -^dr + \Si-^zd,8i + ^dl (5.4) 



it follows from Eq. (3.1) and the assumptions in Lemma 5.1 that 



(1 - /i) - (2Pi(«) + {P2{s) - h)^ \dl^ 



< C{bo,qo){^\V^S2^\ + ^\V^ip\). (5.5) 



Moreover, by the assumptions in the Lemma 5.1 and the expression of /j(l < i < 3), one has 

3 

J2\fi\<C{bo,qo)e. 

i=l 

Combining this with — 2Pi{s)s + ^2(5) < from Lemma 2.2 and (5.4)-(5.5) yields 

r|(9,V,C<^, 92^)1 < C{bo,qom.S<f\ + iVx^l). 

Together with 5*2 = dg and the change of coordinates (1.5), this concludes the proof of Lemma 5.1. □ 

Now we begin to establish the higher-order energy estimates. 

Since the vector fields S are tangential to the boundary r = bgz, drS'^ip = b^d^S'^ip holds on r = b^z in 
view of the boundary condition (3.3). Therefore, one can apply Theorem 4.1 directly to S""<yj (0 < m < fco — 1) 
to obtain: 

Lemma 5.2. Let the assumptions of Theorem 4.1 be fulfilled. If (p is a C''" [Dx) -solution to problem (3.1) 
with (3.3) and (3.6)-(3.8), where ko > 7, then, for < m < ko — 1 and n < —1, 



j f I V^5™(p|2(r, r, e) drde + C2 /// V^S*"'^!^ dzdrde 

J JboT<r<x{T,e) JJJdt 

+ C3 jj^ z''+\d,S'"ipf dS + Cijj z'^+^^idgS'^ipf dS 

< III CS""^ ■ MS'^ip dzdrde + C5 H ((d.S^'ipf 

+ (5,5"V)' + \ideS"'iP)A (1, r, 9) drdO + Ce // z''+\BoS"'^f dS. (5.6) 



Here, the constants C-i, 1 < i < 6, are given in Theorem 4.1. □ 

As in Theorem 4.3, one needs to control the term / / z^^^{BoS"''(p)'^ dS in (5.6) to obtain related higher- 



order weighted energy estimates. In addition, the term JJJj^ S''"'Cip ■ A4S"''ip dzdrdO appearing in the right- 
hand side of (5.6) has also to be payed attention to. This will produce some boundary terms on the shock 
surface and conic surface, respectively, by integration by parts. 

We now establish: 

Theorem 5.3. Let (f e C^°(p^) and £.{z,e) G C''°{[Q,2-k\ x [1,T]) be solutions to (3.1) with (3.3) and 
(3.6)-(3.8), where ko > 7, and 



E ^''\di'dH\+ E ^'|VL+Vl<M£, (5.7) 

o<(i+i2<[-^]+i o<;<[^]+i 

where M > is a constant. Then, for sufficiently small e > and /i < — 1, 

II E T^'+''^'\^i^'^\\T,r,e)dS+ III z^'+f^lVi+'^l^zdrde 

JJb„T<r<>c{T,e) o<(<fc„_i JJJDt o<i<feo-l 

+ // E ^''+^+'|VL+VPrf5<C(6o,9o)£^. (5.8) 
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Proof. Note that on the shock surface F 

SoS-^V = SvBo^ + Bo(5"' - S^)^ + [Bo, S^] 

and ^ ^ 

S^^ = -{f Ur{so+tOdt)zS?'^+[S?\-{f Ur{so + tOdt)z]^; 
Jo Jo 

here and below [•, •] stands for the commutator. Together with Lemma 5.2 and the initial condition (3.8), this 
yields 

T'^+i // \V^S"'<p\'^{T,r,e)drde+ [ [ [ z^'\V^S"'<fi\'^ dzdrdO + [ f z^'+'^\Va^S"'(p\'^ dS 

J JboT<r<x{T,e) JJJdt J Jtt 

4 

<C{bo,qQ)e^ + C{bo,qo)Y,Ik, (5.9) 

i=l 

where 



((( CS"'ip ■ MS"''<pdzdrde, 

J J J Dt 

1 2 
l2 = f[ 2^-' [SF', ( [ Ur{so+tOdt)z]^ dS, 

J Jtt Jo 
h = I! z^'+^lBoiS^ - S^)^\'' dS, 

J JTt 

h = If z^'+^\[Bo,S^]ip\^dS. 

J J Ft 

To prove Theorem 5.3, we treat each Ik in (5.9) separately. With this aim, we divide the rest of the proof into 
four steps. 

Step 1. Estimate of Ii. First we derive an explicit representation of CS"''ip. Due to <S'i(-) = <S'2(-) = 

and Si(i) = -i,52(i)=0, 
r r r 

CSi(p = SiC(p — 2C(p, CS2'p = SiL(p. 

By induction, one obtains 

£S""<^ = S""£<^ + Ci^S^C^. (5.10) 

0<(<m-l 

Hence, it follows from Eq. (3.1) and (5.10) that 

5'"£^ = /ii+/i2+/i3, (5.11) 

where 

hi = his, v.ip)dls"'ip + Ms, v,ip)dls"'^ + /3(s, v.^)d^,s"'^ + ^h{s, v.^)^^^™^ 

+ \h{s, V,0)a2,5"'vj + \h{s, V,.9i)52,5"V, 
/i2 = his, V.0)[5-, dl]^ + /2(s, V,^)[5-, dl\^ + his, V,0)[5'", dl\^ + ^Uis, V,^)[5'", S^]^ 

+ Ims, Vx9i)[5", d%]^ + ^Ms, v,^)[5'», d%]ip, 



0<l<m ^li+l2<l,li>i 



^ 0<;<m Ni+/2<i ^ 

Let us stress that only In contains derivatives of (p of order m + 2, while I12 and /13 contain derivatives of <f 
of order at most m + 1 which are thus lower-order terms. 

^From the expressions of /i, 1 < i < 7, in (3.2), the inductive hypothesis (5.7), and Lemma 5.1, one obtains, 
for m < ko — 1, 

\h2\ + \h3\<C{bo,qo)s J2 IVx^Vl', 
0<i<m 

which implies 

[[[ {Ii2+Ii3)-MS"'(pdzdrd0 <C{bo,qo)s fff z" ^ \V:cS^ip\^ dzdrdO. (5.12) 

JJJDt JJJDt o<l<m 

Next, we treat the troublesome term / / / In ■ M.S'^ip dzdrd6. In view of the formulas dy^vdy^v = 

J J J Dt 

h^y,{{9y,vf), d^^vdy^v = dy^idy^vdy^v) - ^dy^iidy^v)^) and d^^y^vdy^v = ^{dy^idy^vdy^v) + dy^{dy^vdy^v) - 
(^y^i^^yi'^'^yi'^)} 1 one has 

111 • MS'^if = dJl^ + dril + delfi + Iti, 

where 

III = Ifiz'-raid^S'-^f - i/3Z^ra(a,S"V)' " ^hz^raideS"'^f 

- ^Uz^radeS^ipd^ip + hz^+^bd^S^^drS^^p + \f2Z^+H{drS'^<pf 

+ ^fsz^+^bdeS^'ipdrS"'^, 
2r 

I'll = i/2^'^ra(a,5"'v)' +/3^^raa,.S''"(/)9,5''"(/> + ^/6Z^raae5"V9z5''V 
- \fiz^+'b{d,S"'^f + \f3z''+\drS"'^f - l^f^z^+H{deS^^f 

+ ^hz^+^bdrS^'ipd.S"'^, 

2r 

III = -\-f4z''radeS'"pd,S"^p + i/5^^ra(a,5™v')' + hez''rad,S"'pdrS"'p 
r {r 2 2 

+ ^f4zl'+^bdeS"'<pdrS"'ip + ^f5z''+^bd,S"'ipdrS"'<p + ^fez''+^b{drS"'ipfy 

|-ia,(/iz''ra) - ^dr{f2z''ra) - ^de{h'5Z^ra) + is,(/i^^+i6)|(^,5"V)2 

+ i^-drihz^ra) - \de{^z^ra) - d,{hz^+'b) - ^a,(/5^''+^6)|9.5™9i9.5'"^ 
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1^ 2 r 2 r J 

+ (]:d,i-hz^ra) - \deihz^+'b) - ^d,{-hz^'+'b)}drS^^deS"'^ 

Due to assumption (5.7) and in view of the properties of /», 1 < i < 7, one has 

3 

1=1 

which yields 

/// hi-BoS"'<fidzdrde 

J J J Dt 

< C{bo, qo)e( III z'^lV^oS'^'fil^dzdrde + T'^+i / / | V,,5'"(p|2(r, r, 6) drdO 

VJJJDt J J boT<r<x(T,e) 

II \S7,S"'ip\\l,r,9)drde+ II z^^+'\V ,S^^\^ ds] + 1 1 (boll 

J Jbo<r<x(l,e) J JVt J JJBt\ 



11 - l!i 1 dS 



The boundary condition drS"^(p = bodzS"^<p on Bt, together with condition (4.6) and /6|s=6o = ^o/5|s=6o5 
gives by a direct computation 

and further 

jjj hi-BoS"'<fidzdrde 



I Dt 



< C{bo, go)e| /// z^\V^S"'^\''dzdrd0 + T^+'' 1 1 |V^5'Vl'(r, r, O)drd0 

VJJJDt J JboT<r<x{T,e) 

+ 11 \V,S"'ip\\l,r,e)drde+ II z^+'\V,S"'<p\''ds\. 

J Jbo<r<Y(l.e) J Jtt ) 



/feo<r<x(l,0) -I -ITt 

Combining this with (5.11), (5.12), and the initial condition (3.8) yields 



|/i| < C(6o, go)£ \ III ^"iVxSVP dzdrde + T^+i / / \V^S^<fWT, r, 6) drdO 

+ jj^ ^'^+i|V,5V|'rf5 + e2|. (5.13) 
In particular, for m = 0, one has 

/// C^p- MifdzdrdO < C(6o, go)e( /// z''\V dzdrdO 

JJJdt \J J J Dt 

+ T^+' II \V,<f\''{T,r,e)drde+ II z^+'\V,^\'' dS + . 

J JboT<r<x(T,e) J Jtt ) 
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Substituting this into (4.27) yields 



f f |Vx¥'P(T, r, e) drde + C,, /// ^''iVx'^P drdOdz 

J JbaT<r<x(T,e) J J J Dt 

J JFt 

Step 2. Estimate of l2- The vector fields Sr coincide with the vector fields 5 on F and 
computation yields 

■'O mi+m2=m,mi>l ^''O ' 

^ l<J<m ^ 

Combining this with (3.7) gives 

[S^,{f Ur{so+tOdt)z]^ <C{bo,qo)z(Ms V I^Vl + IClY 

C{bo,qo)(Me J2 II z^+^\W.S'<f\^d8+ ff zi^+^^l'ds). 



and further 



^2 < 



Step 3. Estimate of I^. One has that 

BoiS"" - S^)ip = [Bo, ^™ - S^]ip + (5™ - S^)Bo^. 
It follows from assumption (5.7) and a direct computation that 



[Bo,S"'-S^]^ 



<C{bo,qo)Me Y ^'|VL+Vl 



o<;<m-i 



and 



\{S"'-S?)Bo^\<Cibo,qo)Ms ^ ^'|VL+V|- 



0<i<m 



Combining (5.16)-(5.18) yields 



\h\<Cibo,qo)Me Y // z^'+^+'\Vl+'^\'dS. 

Step 4. Estimate of I^. As in Step 2, one has 

[Bo,S^]<fi=[Bo,S"']^ onr. 

Due to 



0<i<m-l 
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and 

+ M2^'C = S'nii, V,^) on T, 

one then arrives at 



and 



o<fe<( o<fc<; 



|/4|<C(6o,go)( E // z^+'\^,S'^\^ds + ff z^^+^eds). (5.20) 
Substituting (5.13)-(5.15), (5.19)-(5.20) into (5.9) and applying Lemma 5.1 yields 

// E T^'^'^'\Vr^\\T,r,e)dS+ [[[ z''^'\yi^'^\'dzdrde 

J JboT<r<x{T,9) o<l<ko-l "'^^ 0<;<feo-l 

+ // E z'^+^+^\Vl+'cp\'dS<C{bo,qo)(ff z^-'\zedS + eA. (5.21) 
-^^T o<(<fe(,_i \JJrT / 

On the other hand, by the boundary condition (3.7), Lemma 4.2, and (5.14), one has 

ff z''-M\^dS<C{bo,qo) ff z^-^^^dS<C{bo,qo)e^. 
Together with (5.21), this yields 

// V T^'+^+^\V'+^^WT,r,e)drdO+ [[[ V z^'+f^lVl+^^l^ dzdrdO 

JJboT<r<xiT,e) o<i<feo-l JJJdt o<l<ko-l 

+ 11 E ^''+''+'|vL+vprf5<c(6o,?o)£^ 

J J'^T 0<i<feo-l 

which shows that the conclusion of Theorem 5.3 holds. □ 

§6. Proof of Theorem 1.1. 

Based on the higher-order energy estimate established in Theorem 5.3, we now prove the global existence 
of a conic shock wave, as asserted in Theorem 1.1, by using a local existence result and continuous induction. 
For any given zq > 0, the solution to (1.8) with the initial data given on z = zq and boimdary conditions 
(1.9)-(1.11) exists in an interval [zq, ^^o + C] for some > by the local existence result of [24] or [11, Appendix] 
provided that the initial data is smooth and satisfies the compatibility conditions. Moreover, if the perturbation 
of the initial data given on z = zq is 0(e), then the lifespan of the solution is at least C/e, with some C > 0. 
Therefore, as long as one can establish that the L°°-norms of (p, ^, and their derivatives decay with a rate in 
z, the solution can be continuously extended to the whole domain. That is, by the local existence result and 
the decay properties of the solution one obtains the uniform boundedness of 99, ^, and their derivatives, and 
then one extends the solution continuously from zq < z < zq + to zq + < z < zq + 2^, with C > being 
independent of zq. Hence, the key to proving Theorem 1.1 is to establish the decay of the L°°-norm of (p, ^, 
and their derivatives. 

It follows from Sobolev's embedding theorem (see also [11, Lemma 14]) and the assumptions of Theorem 5.4 
that, for boz <r< x{z, 0) and 1 < z < T, one has 

0<;<feo-3 ■^•^6o2<r<x(^,0)o<i<feo-l 
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On the other hand, (5.8) shows that 

// J2 l^'^x+Vl' drde < C{bo, qo)e^z-^-\ 

J JboZ<r<x(z,e) o<(<fc„_l 

Hence, ^ |z'V|,+Vl^ < ^(60, go)^^^"''"^ for boz < r < xiz,0) and 1 < z < T. For ko > 7, one 
o<;<feo-2 

has ^ |z'V^+Vl < C(6o,go)£-2"^"^- In addition, due to fco - 3 > [y] + 1, Eqs. (3.6) and (3.7) yield 

i<[^]+i 

z^^\dl^dg^^\ < C(6o5 ?o)£-2~^~^- If one now chooses fi G (—2,-1), then it follows by continuous 

0<il+i2<[^] + l 

induction that the proof of Theorem 1.1 is complete. □ 

Appendix 

Proof of Lemma 3.2. First we estimate ^5(5) and Pi{s). By (3.2), Lemmas 2.1-2.2 and Remark 2.2, then 
1 / 7-1 2^ /„^.. /„^ , 7 + 1. 



+ SU^{s)u',{s) + ^^^Y^M^(s) - ^c2(p(s))^ 

l + 0((6ogo)~^) + 0((6o'Jo)"^)/ 7-1,2 ^o^o Mo 



-bi- 



(T+bJF 2(1+62) {'mo) \ u V uy 

7+1 J, baQQ qo . qo 6o«o ,7-1 {boqo)^ 7-I 



+ ^0 • rTT2 • , ^2^2 + -^TrTT2^^ " TTTTT^tCMo) 



2 ^ l + 6g (1 + 62)2 - 1 + ^,2 (1 + ^,2)2 2 (1 + 6^)2 4{l + biy 

' ■0((6o9o)-') + 0((Mo)"- 



-2^bl{l + bl) / 



Next we derive an expression for P{(s). Recall from (3.2) that 

Uz{s)Uris) 



Piis) 



Then 

u'^{s)Uris) + U;,{s)u'^{s) U;,{s)Ur{s){2u^{s)u'^{s) - ^7(7 - 1) p'> '"^ {s) p' {s)) 



u^,{s)-cHp{s)) iul{s)-c'^{p{sW 
It follows from Lemmas 2.1-2.2, Remark 2.2 and a direct computation that 



p ^ (TfIft^CI + Q((Mo) — )) - (T^T^(1 + 0((Mo)-^)) 

' ' ^_^(l + 0((6o9o)-^))-2(^^(Mo)2(l + 0((Mo)-2) + 0((6o«o)-^)) 

i^TTt (1 + 0((Mo)-') + 0((Mo)-^))(2^ + 0((6o9o)-')0((Mo)-^)0(^) 



+ O((6o?o)-^)) - 2(I^fe) (Mo)2(l + 0((6o9o)-2) + O((6o?o)-^)) 

(-1 + 62) (i + 0{{boqo)-^) + Oiibaqo)-^)) 2bl (l + O((6o9o)-') + O((6o9o)"^) 



(1 - ^-^hlH + bl)){l + 62) (1 - 2^62(1 + 62))2(l + 62) 



,!_2:iiA6 / 

2 "0 2 "O / 1 , ^//. „ x-2- 



1 -I- 2^62 - ini^e 
^ " ^ " 'l + 0((6oQo)-^) + 0((Mo)-^) 



(1- V&J(l + 62))2(i + 62) 
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Other estimates in Lemma 3.2 can be carried out analogously in terms of the expressions for Pi, 1 < i < 5, 
in (3.2) and Lemmas 2.1 and 2.2; we omit the details here. The proof of Lemma 3.2 is complete. □ 

Next, we provide the proof of Lemma 3.3. 

Proof of Lemma 3.3. Prom the expressions for Bi, i = 1,2, 3, in (3.5) and Lemmas 2.1 and 2.2, one has 

( 2A7a+b,i) )^(^ogo)^ bpgo ( (bogo)^ go 6^go 
2(^^%(Mo)2 i + ;,2^ (1 + 62)2 1 + 621 + 62 

= TT6| (^I^(rT6i)) + + o((6o.o)-^)) . 

Similarly, 

= (^^5:J(^)^(Mo)^.o[^(i + o((Mo)-^^ 

= -RrW(^^WT6i))'" + ^o((Mo)-) + ^o((Mo)-^)). 

On the other hand, expressions for fXi(i = 1,2) can be obtained from the estimates of which completes the 
proof of Lemma 3.3. □ 
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